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ABSTRACT 

We re-assess the question of a systematic time delay between the formation of the 
progenitor and its explosi on in a type la super nova (SN la) using the Hubble Higher-z 
Supernova Search sample l)Strolger et al]l2fl04l) . While the previous analysis indicated 
a significant time delay, with a most likely value of 3.4 Gyr, effectively ruling out 
all previously proposed progenitor models, our analysis shows that the time-delay es- 
timate is dominated by systematic errors, in particular due to uncertainties in the 
star-formation history. We find that none of the popular progenitor models under 
consideration can be ruled out with any significant degree of confidence. The inferred 
time delay is mainly determined by the peak in the assumed star-formation history. 
We show that, even with a much larger supernova sample, the time delay distribu- 
tion cannot be reliably reconstructed without better constraints on the star-formation 
history. 

Key words: supernovae: general — cosmology: observational. 



1 INTRODUCTION 

Type la supernovae have been used extensively as standard 
distance indicators and have provided th e best evidence to 
date for an acceleration of the Universe iRiess et al.lll998t 
iPerlmutter et alJll999t iRiess et alJl2004h . Future missions, 
e.g. GAIA and SNAP, will greatly increase the number of 
detected SNe la and significantly reduce the statistical er- 
rors in the determination of cosmological parameters. How- 
ever, the nature of the progenitors of type la supernovae is 
still unknown a nd the empirically calibrated Phillips relation 
JPhillipsl fl993l is not fully understood physically. 

Several progenitor scenarios are under discussion, but 
there is no conse nsus due to uncertainties in the evolu- 
tiona ry ^ processes llHachisu fc Nomotolll99 fil: lHac hisu et all 

;ll 119971 lLanger et alJ 



between the redshift dependence of the supernova rate 
(SN R) and the sta r -forma tion history (SFH). 

IStroleer et alJ l|2004l . hereafter S04, aimed to detect 
this difference by studying the distribution of 25 high-z 
SNe la in the Hu bble Higher-z Supernova Search sample 
l)Riess et alJl2004T) . Their approach was to infer the mean 
time delay of the distribution using a Bayesian analysis, 
which assumed different parametrized time-delay distribu- 
tions and adopted the star-formation history (SFH) from 
iGiavalisco et alJ i2004) . hereafter G04. They concluded that 
mean time delays shorter than ~ 2 Gyr ought to be excluded 
with a 95 per cent confidence level, ruling out essentially all 
progenitor scenarios currently u nder discussion . In a re cent 
re-assessment of the constraints, IStroleer et alJ l|2005fJl ob- 



Li fc van den Heuve' 



lHan fc Podsiadlowskill20oir an d the explosion mechanism 
llHillebrandt fc Niemeverl 1200(1 iRonke fc Hillebrandd Eoost 
iGamezo Khokhlov fc Oranl l2005lh One of the signatures 
of the various scenarios is the distribution of time delays 
between the formation of the progenitor systems and their 
explosion, which could give rise to a significant difference 
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tained a 95 per cent lower limit ranging from 0.2 to 1.6 Gyr 
for different time-delay distributions. In the corrected best- 
fitting model, the 95 per cent confidence interval ranged 
from 1 to 4.4 Gyr with the most likely value at 3.4 Gyr. 

Unlike core collapse SNe (SNe II, Ib/c) that origi- 
nate from massive progenitors with relatively short main- 
sequence (MS) lifetimes (~ 3 — 20 Myr), SNe la are believed 
to be thermonuclear explosions of white dwarf stars (WDs) 
whose progenitors have MS lifetimes ranging from ~ 30 Myr 
to several billion years. This implies a minimum time delay 
for SNe la of the order of ~ 30 Myr. 

Most of the SN la progenitor scenarios that have been 
proposed involve mass transfer on to a CO WD in a binary 
system, either through the expansion and Roche lobe over- 
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flow of an evolved companion (single degenerate [SD] sce- 
narios) or through the slow release of gravitational waves, 
orbital shrinking, Roche lobe overflow and merging of a com- 
pact double WD system (double-degenerate [DD] scenarios) . 
Both scenarios have associated time-delay distributions that 
have been estimated with binary population synthesis codes 
(BPS), where the properties of binary systems are followed 
from their birth up to the explosion stage through the many 
different evolutionary paths. Independently of the particu- 
lar treatment of the binary interactions, the resulting time- 
delay distributions differ considerably since their character- 
istic time-scales have different origins. 

The SD scenario is controlled by the process of 
mass accretion, which has to occur at just the correct 
critical rate in order to allow the growth of the mass 
of the companion WD up to the Chandrasekhar limit 
iNomoto fe Kondcl ll99lT) . The dominant evolutionary 
path seems to occur via the accretion of matter on 
to a CO WD 
CO 



from a slightly e volved MS star, the 



WD 



MS 



1992: iRaDDaDort. Di Stefano. & Smith! 
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In this channel, the ac- 



cretion rate is determined mainly by the mass of the donor 
star, which must lie in a narrow range in order to satisfy the 
required accretion-rate constraints. As a consequence, the 
distribution of MS lifetimes and the time-delay distribution 
of the channel are relatively narrow, peaking at ~ 670 Myr 
and rapidly becoming negligible after ~ 1.5 Gyr. 

Although recent simulations have suggested that other 
evolutionary paths within the SD frame work are of minor 
importance llHan fc Pods iadlowski 120041) . it is quite possi- 
ble that their contribution has been underestimated. This 
is particularly important for the CO WD + RG - SD sce- 
nario, where a red-giant (RG) star acc retes matter on to a 
CO WD star (lHachisu fc Nomotolll996l) : in this channel the 
time-delay distribution extends up to several Gyr. 

The DD scenario Jlben fc Tutukovl 11984 IWebbinkl 

1984) . in contrast, is controlled by the time that it takes 
for the binary system to coalesce, which depends roughly on 
the four th power of the separation o f the double-degenerate 
system dShapiro fc Teukolskvlll983) . As a result, the time- 
delay distribution can be described by a low time-delay cut- 
off (~ 30 — 100 Myr) and an approximately power-law de- 
cline up to the age of the Universe. The lower time-delay 
cutoff can be explained by the time required to form the 
most massive degenerate systems with the shortest MS life- 
times, whereas the power-law tail can be explained by the 
power-law relation between coalescence time and separation 
of the double-degenerate systems. 

However, the expected accretion rates in the DD sce- 
nario are a problem: they are so high that present calcu- 
lations suggests that this leads to accretion-induced col- 
lapse (AIC) and the formation of a compact object rather 
than a thermonuclea r explo si on jNomoto fc Ibenl Il985t 
Saio fc Nomotol Il985l Jl998t iTimmes. Wooslev. fc Taaml 



1994: lNomoto fc Kondo 1991 



Therefore, the currently generally most favoured pro- 
genitor scenario is the SD scenario. Because the seemingly 
dominant evolutionary path of this channel would need to 
be discarded if the mean time-delay were found to be higher 



than 2 Gyr, it is important to confirm the significance of the 
S04 results. 

In this work we have studied the SN la time-delay distri- 
bution using the sample of S04 and the same basic analysis, 
but introducing alternative SFHs found in the literature, 
avoiding binning effects as much as possible and using a 
Goodness of Fit (GoF) test that is generally recommended 
for small samples. We discuss the data and analysis in Sec- 
tion^ show the results and Monte Carlo simulations in Sec- 
tions|H]and|3]and discuss their significance in Sections |5| and 
HJ Throughout his paper, we adopted a value for the Hub- 
ble constant of Hq — 70 km s _1 Mpc -1 , present ratios of 
matter, curvature and dark energy density over the critical 
density of Qm = 0.3, £Ik = and = 0.7, respectively, 
and a 'dark energy' pressure over density ratio ('equation of 
state') of w = — 1. 



2 DATA AND ANALYSIS 

The analysis is bas ed on the Hubble Higher-z Super- 
nova Search sample (iRiess et al.ll2004t iDahlen et al.l l20o3 : 
IStrolger et al.ll2004h . which contains 25 SNe la found in the 
GOODS field (13 in the Hubble Deep Field North, HDFN, 
and 12 in the Chandra Deep Field South, CDFS) in the red- 
shift range 0.21 to 1.55. The SNe were discovered in four dif- 
ference images that were produced by observing both fields 
five times in intervals of approximately one month, compa- 
rable to the typical duration of the main SN la light curve 
peak. 

To infer the underlying time-delay distribution of SNe 
la, S04 compared the observed redshift distribution in the 
sample with a parametrized predicted distribution, derived 
from the G04 SFH convolved with three alternative time- 
delay distributions. 

Each distribution was parametrized by its mean time 
delay, which was recovered using a Bayesian analysis. Among 
these, the distribution that best fit the data was a 'nar- 
row Gaussian', which after being corrected was centred on 
3.4 Gyr with a FWHM of ~ 1.5 Gyr. The 95 per cent con- 
fidence interval for the mean time delay ranged from 1.0 to 
4.4 Gyr. The alternatives 'wide Gaussian' and e-folding dis- 
tributions had a mean time delay above 0.2 and 1.6 Gyr, 
respectively, with more than 95 per cent confidence. 

Only the shapes of the distributions are compared, i.e. 
the analysis is scale-free, and the associated efficiencies of 
SNe per unit formed mass are calculated later by normal- 
ising the models to the SN numbers and are not used to 
constrain the models. This means that the sample must 
ideally span a redshift range that includes both the rising 
and declining parts of the SNR, i.e. where the SNR is not 
appro ximately linear. A recent study (see iBarris fc Tonrvl 
2005, Fig. 6-8) could not fully exploit information on the 
SN redshift distribution because their sample did not reach 
to a sufficiently high redshift (z > 1), as the authors indi- 
cate in the text. A simil ar situation is found in the work of 
iGal-Yam fc Maozl (120041) . 

Thus, because the Hubble Higher-z Supernova Search 
sample is the deepest SN sample available, it is the most 
suitable one for constraining the time-delay distribution of 
SNe la. 

However, the formal errors quoted in S04 reflect only 
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the limited size of the sample and not other systematic un- 
certainties, such as those associated with the SFH. 

In the following Sections l2.1l to l!01 we introduce the for- 
malism that gives the SNR, the number of detected SNe per 
unit redshift and the control times used in the derivations. 
In Section Li.4l we discuss alternative time-delay distributions 
and in Section f2.5l alternative SFHs. The Bayesian analysis 
is described in Section 12.61 and further modifications con- 
cerning binning effects and the GoF test are discussed in 
Section |2~7I 

2.1 The SN la rate 

The rate of SNe la per unit time per unit co-moving volume 
(SNRia.) is given by the star-formation rate per unit time 
per unit co-moving volume (SFR) convolved with the nor- 
malised distribution of explosions per unit time of the pro- 
genitor channel (the time-delay distribution, <f>), and mul- 
tiplied by the number of SNe per unit formed mass (the 
efficiency, v). We assume that neither v nor <f> evolve with 
redshift: 

SNRi^z) = v j SFR(t') (j}(t - t', r) dt' , (1) 

where t = t[z), r is some characteristic time-scale defined 
in Section 12.41 and Zr is the redshift associated with the 
time when the first stars formed, approximately the epoch 
of reionisation. We assumed zr = 10, as in S04. 

2.2 Distribution of detected SNe 

The number of detected SNe la per unit redshift interval 
(nia) is given by the multiplication of the rate of SNe la 
per unit time per unit co-moving volume (SNRia), a time 
dilation factor, (l + z)" 1 , the control time of the survey (t c ) 

and the volume per unit redshift being surveyed, r— r Acj: 



Time delay distributions 



dzdu 



. SNR^z) dV{z) . 

ni a («) = — — — — t c (z) -^f 1 Au), 



1 + z 



dzdui 



(2) 



where in our cosmology the volume derivative formula sim- 
plifies to: 



Ic—'t—^-, where 



-1/2 



(3) 
(4) 

(5) 



dV 

dzdtu dz 

d c = cHg 1 f du [(1 + u) 3 n M + fi A ] " 
Jo 

and hence, 

= cff-i [(i + z fn M + n A ] ~ 1/2 d 2 c (z 

2.3 The control time 



The control time can be understood as the total observ- 
ing time multiplied by the probability of detecting a SN 
at a given redshift. We used the same values as S04, that 
were calculated taking into account the expected extinction, 
spectra, light curve shapes and peak magnitude dispersion of 
SNe la, the way each field was revisited, and the efficiency of 
the detection algorithm (but see Section ^. 7. H , The control 
times are defined by 



SD scenario 
DD scenario 



— Narrow Gaussian (1 Gyr) 

— Narrow Gaussian (4 Gyr) 
e-folding (1 Gyr) 
e-folding (4 Gyr) 




Time [Gyr] 



Fie 



1. 



Theoretica l time delay distributions 
iHan Pod siadlowskil l2004h compared to parametrized time 
delay distributions used in the analysis. The best-fitting model 
in S04 corresponds to the 'narrow Gaussian' distribution with a 
mean time delay of 4 Gyr 



Time [Gyr] 

10 9876543 2 1 




0.5 1 5 10 

Redshift 

Figure 2. Histogram of progenitor formation redshifts for the 25 
SNe in the Hubble Higher-z Supernova Search sample assuming 
unique time delays. Note that some SN progenitors must have 
originated at extremely high redshifts if the time delays were al- 
ways high. 



t c (z) = 



P(t\M x ,A x ,z) P(M X ) P(A X ) dM x dA x dt, 



(6) 



where P(t\M x , A x , z) is the probability of detecting a new 
SN at time t, given its rest-frame luminosity and its host 
galaxy extinction and redshift. It depends on the assumed 
spectra through K-corrections, the sensitivity of the survey 
and the efficiency of the detection algorithm. 

P(M X ) is the probability of having a given SN rest- 
frame luminosity. It was estimated based on the character- 
istic relation between peak luminosity and light curve shape 
of SNe, and the observed dispersion of SN la peak luminosi- 
ties. 

P(A X ) is the probability of having a given host galaxy 
extinction at the given rest-frame wavelength. It was as- 
sumed to be proportional to e~ Ax . 

For more details see the original discussion in S04. 

2.4 Time delay distributions 

All the time delay distributions were parametrized by their 
mean time delays, r. S04 used an exponential distribution 



© 0000 RAS, MNRAS 000, 000-000 



4 F. Forster, C. Wolf, Ph. Podsiadlowski and Z. Han 



and two families of Gaussian distributions whose width scale 
with the mean time delay. The e-folding distributions are 
given by: 



-t/7 



(7) 



The two alternative Gaussian distributions are grouped into 
the families of 'narrow' (ov = 0.2r) and 'wide' (ov = 0.5t) 
distributions, of the form: 



1 



(8) 



The previous distributions are defined only for positive val- 
ues. However, negative time delays must be allowed in order 
to avoid statistical bias in a small sample and to get confi- 
dence intervals that do not artificially discard short time de- 
lays. Moreover, a preference for negative time delays would 
signal SFHs that peak too late in time. Thus, we consid- 
ered a fourth time delay distribution, a Gaussian distribu- 
tion with fixed width (a = 0.5 Gyr) that allows either for 
positive or negative time delays: 



cj>(t,T) = 



(9) 



We also added a log-normal distribution, which is associated 
with processes where the source of uncertainty has multi- 
plicative effects rather than additive ones, as is the case for 
Gaussian distributions. The best-fitting models to the the- 
oretical time delays were in most cases log-normal distribu- 
tions, whose width a, in units of log(t), was kept fixed and 
determined by the best-fitting model of the theoretical time 
delays: 



" 2 ln(<r) a Z 



(10) 



The theoretical time delay distributions from 
lHan fc Podsiadlowsk i (2004) were also examined with 
a GoF test. They were produced assuming either the CO 
WD + MS - SD scenario or the DD scenario with different 
binary evolution parameters. In Fig. the theoretical time 
delay distributions of the SD and DD scenario together 
with two time delay distributions with different mean time 
delays are plotted, including the best-fitting model from 
S04. 



2.5 The Star-Formation History (SFH) 

Because we consider alternative prescriptions for the SFH 
that have incomplete redshift information, further compli- 
cations arise. The ideal redshift coverage of the SFH should 
range from zero to the redshift of the first star formation, 
farther than the highest-redshift object presently known in 
the Universe. This is a consequence of the redshift range 
of the detected SNe and the long time delays that have to 
be considered. If high-redshift SNe only exploded after long 
time delays, their progenitors would need to form at red- 
shifts up to ~ 30, as Fig. |5|shows. 



If the determination of the SFH did not cover the re- 
quired redshift range, we used as an approximation either a 
power-law extrapolation in time or a scaled version of the 
optical-UV derivation. We found that the method is not 
very sensitive to this approximation when the position of 
the peak of the SFH is well constrained, since it is the dif- 
ference between the peaks of the SFH and the SNR that 
mainly constrains the best-fitting models. The alternative 
prescriptions of the SFH we have used are the following: 

• The SFH bv lGiavalisco et alJ i2004\ . G04. We have used 
continuous approximations for the extinction corrected and 
not corrected models, inferred from deep optical-UV obser- 
vations of galaxies in the GOODS field. Both versions differ 
by a factor of ~ 3 in the redshift range of interest. The 
continuous approximations are the ones used in S04, which 
peak at z ~ 2.7 in the extinction corrected model (Ml) and 
at z ~ 1.8 in the model that is not corrected for extinction 
(M2). 

• The best-fitting model of ICharv fc Elbazl fcOOlft . here- 
after CE01. It was derived from the integrated cosmic in- 
frared background (CIRB) and covers the redshift range 
from to 4.5. At z > 4.5 we tried a power-law extrapola- 
tion in time or a scaled version of G04 at z > 3. Mainly 
because the peak of the SFH occurs very late in time, at 
z ~ 0.8, we found that the inferred time delays are not very 
sensitive to this approximation. However, a constant SFH 
model is within th e error bars at high-z . 

• The SFH from lHeavens et al.l (I2004F) . hereafter H04, in- 
ferred from the 'fossil record' of stellar populations in the 
Sloan Digital Sky Survey (SDSS). We interpolated a Spline 
function to the binned SFH, which peaks at z ~ 0.4, in order 
to obtain a smooth approximation. It is not usually recom- 
mended to approximate data in this way, but we think that 
our approximation preserves the general differences between 
this SFH and the alternative prescriptions. We also tried a 
scaled version of G04 at z > 3, or a constant star-formation 
history for all times, since this SFH is very flat at z > 1 and 
peaks very late in time with respect to the detected SNe. 
Both alternatives gave very similar results to the fit to the 
original binned data. 

Additionally, we considered one of the most recent de- 
terminations of the SFH us ing Spitzer data, presented by 
iPerez-Gonzalez et aljl2005T) . However, both its limited red- 
shift coverage and its dependence on the assumed galaxy 
luminosity function makes the high-redshift extrapolation 
ambiguous. For this reason, we did not try a continuous ap- 
proximation of this SFH in the analysis, although it must be 
considered a reliable result. The variance between its differ- 
ent versions only demonstrates the persisting uncertainties 
in our knowledge of the SFH. 

In Fig. [jj] we show the four continuous approximations 
of the SFH and the binned SFH from PG05. The contin- 
uous approximations are based on the extinction corrected 
and not corrected SFH from G04, the best-fitting model 
from CE01 with a power-law extrapolation in time at high 
redshift and a continuous approximation of the best-fitting 
model from H04, which is constant at high redshift. It is 
apparent that there is a range of SFHs in the literature that 
do not agree and, importantly, peak at very different times. 
Thus, it is important to understand the systematic errors 
associated with this uncertainty. For a recent estimation of 
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Figure 3. The alternative SFHs that better reflect the diversity of resul ts obtained in the ana lysis. For G04 jGiavalisco et alj |20041. 
Ml is the extinction corrected model and M2 is not corrected. For CE01 ICharv fc Elbaz 2001J), we used a powe r-law extrapolation i n 
time for z > 4.5, but also a scaled version of G04 (not shown here), which gave very similar results. For H04 tHeavens et al.lbootfl , 
horizontal error bars represent bin sizes and vertical error bars (very sma ll), bootstrap errors. In our calculations we chose to represent it 
by a smooth Spline interpolation. The PG05 iPerez-Gonzalez et all2005h SFHs are determined combining optical-UV data with Spitzer 
observations in the GOODS field, assuming different galaxy luminosity functions. The SNe in the Hubble Higher-z Supernova Search 
sample are in the redshift range from 2 = 0. 21 to 2 = 1. 55. 



the uncertainties on the SFH see also Fig. 2 and 4 from 
iHopkins fc Beacoml d200rj| . 



2.6 The Bayesian probability 

Using Bayes theorem with a uniform prior, the probability of 
a mean time delay r with a time-delay distribution <f>(t,T) 
and a SFH in the form of SFR(t'), given the set of SN 
redshifts, {zi}, is proportional to: 

P(SFR(t'),^t,r),r\{ Zl }) cc P({ Zi }\SFR(t'),<t>(t,T),T). 

(11) 

Thus, it is proportional to the probability of the particular 
SN redshift distribution: 

25 

P({ Zi }\SFR(t'), 4>{t, t), t) oc [J ni ^ z ' r )' ( 12 > 

where m a (z, r) has been normalised for every r because the 
analysis is scale free. It depends on r through the time-delay 
distributions of Section [2.41 and equations and [5] 

Hence, for a given combination of SFH, time-delay dis- 
tribution and mean time delay, the predicted number of SNe 
per unit redshift can be expressed as a probability distribu- 
tion in redshift. Subsequently, the probability of the set of 
SNe can be calculated for every r. 



2.7 New analysis 

The main differences between the S04 calculations and this 
work are a result of considering a range of alternative SFHs. 
Further differences are as follows: 



2.7.1 Redshift binning effects 

The main advantage of using the observed redshift dis- 
tribution of SNe la in stead of the corresponding SNR 
iGal-Yam fc Maozl 12004,1 is that the analysis can be done 
in a way that avoids binning and the subsequent loss of in- 
formation. 

Moreover, we have found that the analysis is very sensi- 
tive to the way the volume derivative is calculated in equa- 
tion [5] Because the SN sample is small, binning the data is 
not recommended, and all the calculations should be done 
continuously. Binning can introduce a relative overestima- 
tion of the volume derivative at low redshift, effectively 
pushing the most probable time delays to higher values. As a 
result, the lower limit of the Bayesian analysis can be overes- 
timated by more than ~ 1 Gyr and the peak of the Bayesian 
probability distribution by ~ 500 Myr, according to our cal- 
culations. This is con sistent with the corrected results in 
IStrokeret all J2005rj) . 

However, in order to get a continuous version of the 
control times we have used an interpolation of the values 
calculated in S04 at redshift intervals of 0.2. Recalculating 
the control times with more redshift resolution would be a 
better approach, but we have not tried it in this analysis. 



2.7.2 Goodness of fit test (GoF) 

A maximum likelihood analysis must be accompanied by 
a GoF test to check that the parametrized model has an 
appropriate form to start with, and only then can the 
confidence intervals be trusted at all. Accordingly, a \ 2 
test was used to check consistency between the predicted 
and observed redshift histograms of SNe la in S04. How- 
ever, the x 2 t es t is n °t reliable when the number of ele- 
ments per bin is not greater than five in 80 per cent of 
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Figure 4. Bayesian probabilities and KS test non-rejection probabilities for the case of the 'fixed width Gaussian' and e-folding time 
delay distributions. The probabilities associated with the 'narrow Gaussian' distribution are very similar to those associated with the 
'fixed width Gaussian' distribution in the positive time-delay region. Note the relation between the most probable and best-fitting values 
with the peak of the SFHs in Fig. [3] i.e. later peaked SFHs give shorter time delays. The apparent preference for negative time delays 
in the case of the SFH from H04 may be an indication that this SFH peaks too late in time. See Section|2]for more details. 



the bins ijWall fc Jenkins! I2003T ) . Instead, we have used a 
Kolmogorov-Smirnov test (KS test) as our goodness of fit 
test, which is the recommended test to use when the sample 
size is small and because the analysis is done continuously. 
Furthermore, selecting the best-fitting models with the KS 
test can be used as an alternative parameter estimator. 

2.7.3 Confidence intervals 

In S04 the confidence intervals were obtained starting from 
the mode (the maximum) of the Bayesian probabilities, 
partly because in the original results the probabilities at 
low time delays were negligible. However, because in some 
cases the probability distributions are relatively flat, or not 
negligible at zero time delay, the definition of the confidence 
intervals becomes important. In our calculations, taking 95 
per cent confidence intervals around the maximum vs. from 
the median can make a difference of typically ~ 500 Myr. 

One way to avoid this problem it to use the parameter 
region that is not rejected by the GoF test with a certain 
confidence level. In this approach, we obtain 95 and 68 per 
cent confidence intervals that are unambiguously defined. 

2.7.4 Photometric redshifts 

Because the spectra of SNe are characterised by many 
blended lines broadened by high velocity dispersion, SN red- 
shifts are determined from their host galaxies. Of the 25 
SNe la in the sample, six have only photometric redshifts, 
three of them in each field. We found that the photomet- 
ric redshift of SN 2003al, 0.91 ± 0.2, h as a better estim ate 
in the public COMBO-17 catalogue dWolf et al.ll2004h of 
0.82 ±0.04. Additionally, in lStrolger et alJ J2005alf . the pho- 
tometric redshift of SN 20031u, 0.11±o;n, has a better esti- 
mate of 0.14 ±0.01. 



3 RESULTS 

3.1 Varying the SFH 

If alternative SFHs are allowed, the Bayesian probabilities 
associated with a given time-delay distribution have a wide 
range of preferred time delays. The Bayesian probabilities 
and KS test associated rejection probabilities show a pref- 
erence for values ranging from very long (~ 4 Gyr) to very 
short, and even negative (~ —3 Gyr) if the SFH peaks very 
late in time (see Fig. ■ The relation between the peak of 
the SFH and the peak of the Bayesian probability distribu- 
tion is, to zeroth order, such that later peaked SFHs give 
shorter time delays. 

Interestingly, inspection of Fig. [I] (upper-left panel) 
shows two types of maxima in the Bayesian probabilities: 
one whose position decreases with later peaked SFHs and 
another that is fixed at approximately ~ 3.5 Gyr, even for 
different SFHs. The first peak approximately reflects the 
time difference between the peaks of the SFH and the SNR. 
The second peak reflects the relative absence of SNe at high- 
z. Because no SNe were detected between the epoch of reion- 
isation and z ~ 1.5, or between t ~ 0.5 Gyr and t ~ 4 Gyr, 
the Bayesian analysis marginally favours models that do not 
produce SNe in the first ~ 3.5 Gyr after the assumed epoch 
of reionisation. The upper plot of Figgis illustrative of this 
effect. 

With the current data, it is the first peak which is sta- 
tistically dominant for all the SFHs, but this may change 
with deeper and wider SN surveys in the future. 



3.2 Kolmogorov-Smirnov test 

With the KS test we find best-fitting mean time delays 
and confidence intervals that are free from the problems ex- 
plained in Section [2.7.31 The rejection probabilities for the 
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'fixed width Gaussian' and e-folding time-delay distributions 
are shown in Fig. 0] 

We found that all the combinations of SFH and time- 
delay distributions had an associated parameter region that 
is accepted by the KS test, which validates the use of the 
Bayesian analysis. Additionally, the parameter estimation 
seems robust in the sense that it gives results that are con- 
sistent with what the Bayesian analysis shows. Moreover, 
the addition of this test shows that the negative time-delay 
peak for the H04 SFH is favoured over the long time-delay 
peak (see Fig. , which may be an indication that this SFH 
is not compatible with the SN data. 



3.3 Confidence intervals 

The confidence intervals were defined as the parameter re- 
gion that cannot be rejected with a certain confidence level 
based on the KS t est. We found that only the extinction cor- 
rected SFH from iGiavalisco et alJ ll2004) has a 95 per cent 
confidence lower limit greater than zero, i.e. around 1 Gyr. 
All the alternative SFHs did not result in a lower limit for 
the time delays greater than zero. A summary of the con- 
fidence intervals obtained with the Gaussian and e-folding 
time delay distributions is shown in Fig. 



3.4 Varying the time-delay distribution 

The non-rejection regions of the five time-delay distribu- 
tions tested in this work can be grouped into three fami- 
lies of results: one family associated with the 'fixed width 
Gaussian' (that allows for negative time delays) and 'nar- 
row Gaussian' distributions, another with the 'wide Gaus- 
sian' and log-normal distributions and one associated with 
the e-folding distribution. 

As a general rule, the narrower the test time-delay dis- 
tribution, the narrower the associated Bayesian probabili- 
ties. However, it is in the long time-delay region where the 
changes are more noticeable, as can be seen in Fig.Qj] This is 
because the abrupt transition that occurs in the SFH at the 
epoch of reionisation is reflected in a less smooth SNR when 
narrower time-delay distributions are assumed. Hence, the 
wider the time-delay distribution, the less pronounced the 
second peak in the Bayesian probabilities (see Section 13. II 
and the longer the time delays allowed. 



3.5 Theoretical time-delay distributions 

We performed KS tests of the theoretical time-delay distri- 
bution varying the BPS parameters and the assumed SFHs. 
As a result, if the extinction corrected SFH from G04 is as- 
sumed, the CO WD + MS - SD scenario alone has a 3 per 
cent probability of not being rejected and the best-fitting 
models are obtained for a DD scenario which has a very high 
mass transfer efficiency. Conversely, if any of the theoretical 
models are assumed to be true, the best-fitting models are, 
in almost all the comb inations, obtained with the SFH from 
ICharv fc Elbazl (1200 ll) . In Table0 we show the non-rejection 
probabilities for the different combinations of SFH, theoret- 
ical scenario and BPS parameters tested in this work. The 
BPS parameters are: c yce, the common-envelope ej ection ef- 
ficiency efficiency as in lHan fc Podsiadlowskil i2004h : qrlof, 



1 Bayesian 95-68% c.i., mode, mean-median 
KS test 95-68% c.i., best-fitting model 




Wide Gaussian 



G04(M1) 
G04(M2) 
CE01 
H04 

G04(M1) 
G04IM2) 
CE01 
H04 



LogNormal 



G04(M1) 
G04(M2) 
CE01 Fj 
H04 



Narrow Gaussian 



Fixed width Gaussian 




-2 2 

Time delay [Gyr] 

Figure 5. Summary of the parameter estimation analysis us- 
ing different SFHs and alternative time-delay distributions. The 
Bayesian probability 95 and 68 per cent confidence intervals are 
plotted with the position of the mode and the mean-median in- 
terval, in addition to the KS test rejection probabilities 95 and 68 
per cent confidence intervals and the position of the best-fitting 
time delay. Note that later peaked SFHs show a preference for 
shorter time delays and that wider time-delay distributions al- 
low for longer time delays. H04 gives consistently poorer fits for 
positive time delays than for negative time delays as the use of 
the 'fixed width Gaussian' distribution shows (see Fig.|3]and the 
discussion in Section l2.4l for the use of negative time delays). 



the Roche lobe overflow mass transfer efficiency, and Z, the 
metallicity. 



3.6 Constraints on the SFH assuming the 
theoretical models 

In addition to the previous constraints on the time delay 
distrib ution we followed the approach of iGal-Yam fc Maozl 
ll2004h to constrain the SFH. The method assumes a par- 
ticular time delay distribution and a SFH that consists of a 
broken power-law smoothly joined at the transition redshift 
zq, proportional to ~ (1 + z) a at high z and to ~ (1 + z) 13 
at low z, i.e.: 
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Table 1. Summary of the KS non-rejection probabilities (per 
cent) for different combinations of SFH, theoretical time-delay 
distribution and BPS parameters. Unless stated otherwise, the 
standard parameters are c?ce = 1-0 , orlop = 1-0 an d Z = 0.02 . 



SD scenario - hce=0.75, hrlof = 0.5, Z= 0.02 
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Figure 6. Constraints on the SFH assuming the theoretical time 
delay distributions. We show the KS test marginalised probabili- 
ties of not-rejecting a model with parameters zq and /3 (see equa- 
tion ^^J. The probabilities are marginalised over a in the range 
-2 to 0. We have also marginalised over /3 or zq and found that 
for the purpose of this work, the most important parameter is zq, 
followed by (3 and a in order of importance. 



SFR(z) oc 



+ 



Zq 



5/3 



-1/5 



(13) 



For more details see lGal-Yam fc Maod <l2004l) . We assumed 
the theoretical time delay distributions and performed a 
KS test for different combinations of power-law indices and 
transition redshift, Zq. We marginalised the probabilities on 
a because it was found that the dependency on this pa- 
rameter is weak, which is expected because the SN sam- 
ple contains very few objects at high redshift. The resulting 
probability distributions for different values of zo and /?, as- 
suming either the SD or DD scenario is shown in Fig. 
Virtually all measurements of the SFH suggest SFR(z) oc 
(1 + z) 13 , with (i in the ran ge ~ 2 to 4 iWilson et alJl2002t 
iPerez-Gonzalez et all20o3) . For (3 > 2 we find that the peak 
of the SFH is most likely between z ~ 0.7 and 1.2. However, 
any rejection at 90 per cent confidence is only possible for 
a peak location at z > 2. The flatter the recent decline in 
the SFH (lower /3), the less clearly constrained is the peak 
of the SFH. 
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Figure 7. Systematic errors associated with the choice of SFH. 
We performed Monte Carlo simulations assuming two different 
SFHs: the labelled SFH to derive the mock sample of SNe, but 
the G04 (Ml) SFH to derive the Bayesian probabilities. We show 
the Bayesian probabilities derived from 10,000 mock SNe for one 
combination of BPS parameters in the SD and DD scenarios, 
respectively. The associated systematic error are of the order of 
~ 2 Gyr. 

SD vs DD scenario 
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Figure 8. Systematic errors associated with the choice of BPS 
parameters. The G04 (Ml) SFH is assumed and the BPS parame- 
ters are changed. Because the time delay distribution are not fully 
sampled, the mean time delays are underestimated in some cases, 
but by less than 0.5 Gyr. If some BPS parameters are changed, 
differences of up to ~ 3 — 4 Gyr are found for the mean time 
delay of the DD scenario, while for the SD scenario they are of 
the order of ~ 0.5 Gyr. 



4 MONTE CARLO SIMULATIONS 

In this section we examine whether our approach is subject 
to systematic biases due to the method itself and the choice 
of the SFH. 

To assess the robustness of the method, we performed 
Monte Carlo simulations for each progenitor scenario sep- 
arately, where we took large sets of simulated SNe (typi- 
cally 10,000) to minimize statistical errors. Because the DD 
scenario can extend to times comparable to the age of the 
universe, its time delay distribution cannot be sampled com- 
pletely; we therefore expect that the recovered values will al- 
ways be biased towards shorter mean time delays. We found 
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Figure 9. Redshift histogram of a complete sample of galaxies in 
the GOODS-Southern field from the VVDS survey (gray), all the 
SNe with spectroscopic redshift in the field (white) and only the 
SNe la with spectroscopic redshift (light gray). Bins are chosen 
to show constant co-moving volume. The dashed line shows the 
average galaxy density (assuming no evolution). Several peaks 
are due to large-scale structure. Note that at higher redshifts the 
width of structures appears artificially enhanced by the distorted 
redshift scale. 



the following: (1) such an expected bias is indeed present for 
the DD scenario, but is always below 0.5 Gyr; (2) for large 
sample sizes, the smallest statistical errors are obtained us- 
ing the e-folding time delay distribution, independent of the 
theoretical scenario assumed. 

To quantitatively estimate the systematic errors associ- 
ated with the choice of SFH, we proceeded in the following 
way: (1) in order to minimize the statistical errors and study 
the systematics, we performed Monte Carlo simulations with 
10,000 mock SNe drawn from the theoretical scenarios, us- 
ing the e-folding time-delay distribution in the analysis; (2) 
we produce a mock sample of SNe with a SFH that differs 
from the one used in the Bayesian analysis; (3) a compari- 
son between the different recovered mean time delays then 
gives an estimate of the systematic error due to the choice 
of SFH (see Fig. Q and Fig.|SJ. The result is consistent with 
what can be concluded from Section [3] i.e. that the bias on 
the mean time delay can be of the order of ~ 2 Gyr, even 
ignoring the H04 SFH. 

We then tested the robustness of the algorithm with 
small samples. In order to do this, we performed Monte 
Carlo simulations with 10,000 different sets of 25 mock 
SNe drawn from the theoretical models. We repeated the 
Bayesian and GoF analysis and found the following: (1) 
when the theoretical time delays are relatively short and the 
SFH peaks early in time, or when the time delays are long 
and the SFH peaks late in time, the Bayesian probabilities 
tend to have very positive or very negative skewness and, 
consequently, the mode of the Bayesian probabilities either 
overestimates or underestimates the theoretical mean time 
delays; (2) the distribution of non-rejection probabilities us- 
ing the KS test is fiat, meaning that there is no significant 
bias towards short or long time delays. Hence, we conclude 
that the mode of the Bayesian probability is not a good esti- 
mator for the mean time delays, while the KS test confidence 
intervals are robust even with s mall samples, genera l results 
consistent with the discussion in lPress et al.l lll992l) . 



5 DISCUSSION 

5.1 Large scale structure 

Large-scale structure (LSS) effects are important in small 
pencil beam surveys when studying the time delay of SNe 
la, even when the SFH and the SNR have been measured in 
the same field. Usually, the SFR is measured as an 'instan- 
taneous' observed rate at a particular location in space and 
time, where for this study a prediction of the SNR will be 
needed for comparison. However, this prediction ought to be 
based on the SFR in the same position in space, not only at 
the time when the SNe are seen to explode, but at earlier 
times as well. 

Since the latter figure can not be measured directly, we 
have to obtain the SFR for earlier times by looking at a 
higher redshift, i.e. at a different location in space. Cosmic 
variance and the dependence of star formation on local envi- 
ronment will lead to a SFR measurement that is somewhat 
different from the past SFR of the target location in the 
same field. This SFR variance will appear whether the SFR 
is measured from the same fields or not, since the locations 
in space where the SFR and SNR are measured will be un- 
related. A 1 Gyr time delay already translates into proper 
distance differences of ~ 500 Mpc at z — 1. 

So, if we do not know the individual star-formation his- 
tories of galaxies in the supernova field, we are best ad- 
vised to use our best knowledge of the cosmic SFH, while 
LSS still leaves an imprint on the observed SNR history. 
We investigate the large scale structure in the CDFS where 
the spectroscopic redshif t survey VIMOS VLT Deep Sur- 
vey <lLe Fevre et alJl2004l VVDS) overlaps almost precisely 
with the GOODS field and is complete to My < —19.5 at 
z < 1. In Fig. |5Jwe show a histogram of galaxy redshifts 
with My < —19.5 in bin sizes chosen to contain a constant 
co- moving volume of (25 Mpc) 3 . A non-evolving and homo- 
geneous galaxy distribution should appear flat in this repre- 
sentation. Well-known over-densities or wall-like structures 
are clearly apparent, especially in the redshift range from 0.5 
to 0.75. Two wall-like structures near z w 0.67 and z w 0.74 
are conspicuous in the distribution and ha ve been observed 
in both x-ray an d optical surveys (see e.g. iGilli et al]l2003l : 
IWolf et alJl2004l . for independent confirmation). 

Also shown is a histogram of SNe la with secure (spec- 
troscopic) redshift determinations. SNe with photo-z's only 
have been omitted due to the large error in redshift which 
makes any association with structures in the galaxy distri- 
bution difficult. 

The average bin in the redshift range from to 1 con- 
tains 7.8 galaxies, whereas the average SNe-weighted bin 
contains 16.7 galaxies (2.14 times the normal average), a 
possible indication that the inferred SNR is affected by 
LSS in this field and redshift range. In order to under- 
stand how these numbers depend on the choice of our 
bins, we have repeated the calculations for bin widths of 
(AV) 1/a = (20, 21. .25) Mpc and found ratios of 2.52, 2.38, 
2.15, 2.22, 2.00 and 2.14, respectively. Interestingly, the ra- 
tios that result using only SNe la-weighted bins are even 
higher: 2.87, 2.68, 2.88, 3.29, 2.55 and 2.97. 

The imprint of LSS in the SNR distribution could be 
corrected in the analysis of time delays. If the underly- 
ing mass density on a given direction h were of the form 
p(n, z) — pq(z)(1 + S(h, z)), the over-densities could be cor- 
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rected by multiplying the control times by the same factor, 
1 + S(h, z). However, we do not correct for density variations, 
because we lack a good determination of the LSS with con- 
sistent quality in both fields. 

5.2 SN rates and efficiencies 

In the process of matching the observed number of SNe with 
a model, it is possible to explicitly calculate the SN produc- 
tion efficiencies and the supernova rate (SNR). Not only 
does the shape of the SNR history constrain the progenitor 
models via the delay times, but also the SN efficiency must 
be compatible with the model. 

As a first consistency test, we compare the directly mea- 
sured SNRs with the parametrized estimates. In Fig. 1101 we 
show two versions using the results of the parameter esti- 
mation with the G04 (Ml) and CE01 SFHs. Clearly, the ex- 
tinction corrected SFH from G04 with the 'narrow Gaussian' 
distribution is consistent with long time delays (~ 3.5 Gyr), 
whereas the SFH from CE01 matches best with an e-folding 
distribution of shorter time delays (~ 1.5 Gyr). The second 
alternative is only marginally favoured by the KS test, so 
both possibilities seem equally plausi ble. 

It is important to mention that in lStrolger et al.l <l2005al) 
a deeper SN search in the smaller Ultra Deep Field and its 
parallel fields (UDF/P) resulted in the detection of four ad- 
ditional SNe with z < 1.4. The lack of high-redshift SNe 
is one of the predictions of the S04 best-fitting model and 
hence it supports the result. However, the authors also con- 
sidered the SFH from CE01, but without any time delay, 
and concluded that it is not possible to rule out this SFH 
with more than 50 per cent confidence. 

The SN efficiencies required to explain the observed 
SNR pose a problem for the theoretical SD scenario, which 
produces too few SNe. This is true for all SFHs and combi- 
nation of parameters, and amounts to a shortfall of 4x to 
10 x depending on the assumed SFH. The DD scenario can 
reproduce the required efficiencies for some combinations of 
BPS parameters. However, changes in the little constrained 
binary fraction and mass ratio distributions, or possibly in 
the initial mass function (IMF), may solve this problem in 
the future. 

One additional challenge is that the time-delay distri- 
bution or the supernova efficiencies may evolve with time: 
if the delay distribution is made up of several components 
from different SN production channels, their relative con- 
tribution could change with time; the supernova efficiencies 
themselves may evolve with time, especially considering that 
the accretion processes cou ld depend on environm ental fac- 
tors such as the metallicity llKobavashi et al Jl99dl . The rel- 
ative absence of SNe at high redshift could be a result of a 
metallicity effect, rather than a reflection of the time delays. 
If this is truly the case, our error bars on the time delays 
would be largely underestimated. The use of compatible star 
formation and chemical enrichment histories will be required 
to tackle this problem. 

5.3 Spitzer SFH 

We have shown that it is crucial to determine the SFH that 
is to be used in the analysis. The recent determination of the 



SFH from lPerez-Gonzalez et all (l2005t) , using three different 
extrapolations of the galaxy luminosity function, was shown 
in Fig. This SFH was obtained by combining infrared 
Spitzer observations with optical-UV data of galaxies in the 
GOODS field. The resemblance between this SFH and the 
CE01 best-fitting model justifies the use of the latter. 

Interestingly, the recent near-infrared and 
sub-millimetric determinatio n of the SFH from 
IWang. Cowie fc Bargerl i2005T) peaks at z ~ 1, simi- 
larly to what is obtained in CE01. These new results 
emphasize the persistent uncertainties in our knowledge of 
the SFH and the problem of extinction corrections in the 
optical. 



6 CONCLUSIONS 

We have found that systematic errors associated with the use 
of alternative star-formation histories are co mparable if not 
larger than the statistical errors reported in lStrolger et all 
( 2004). The position of the peak of the SFH was found to be 
the crucial parameter for the recovered time delays: Later 
peaked SFHs result in lower time delays and vice versa. Fur- 
thermore, the confidence intervals for the time delays depend 
on the functional form of the delay distributions assumed in 
the analysis. The use of wider time-delay distributions, in 
particular the e-folding model, gives considerably longer up- 
per limits for the time delays. 

For the data set under investigation, we found that the 
KS test is better suited to obtain confidence intervals than 
a Bayesian analysis. The KS test confidence intervals are 
unambiguously defined, and we have confirmed their validity 
using Monte Carlo simulations. In contrast, the skewness of 
the Bayesian probabilities and the small sample size result 
in somewhat arbitrary confidence intervals. 

A KS test using the shape of the theoretical time-delay 
distributions shows that the extinction corrected model from 
G04 is incompatible with the CO WD + MS - SD scenario, 
although other SFHs are compatible. The DD scenario can- 
not be rejected at 95 per cent confidence with any com- 
bination of SFH and time-delay distribution. Starting from 
theoretical time-delay distributions, they consistently favour 
the SFH from CE01. 

If we wish to constrain the time-delay distribution and 
possibly discard progenitor scenarios from the redshift dis- 
tribution of supernovae, it is of foremost importance to de- 
termine the cosmic star-formation history more accurately. 
Otherwise, the uncertainty in the SFH will continue to limit 
the interpretation of SN data sets of any conceivable size. 

Secondly, we would need a better understanding of our 
progenitor models: if we ignore factors affecting efficiencies, 
their cosmic evolution will make the delay time distributions 
evolve and produce a situation where it is hard to disen- 
tangle these different effects without qualitatively different 
observations. 

Only after having solved these issues would deeper and 
wider surveys help to constrain the delay times and progen- 
itor models by decreasing statistical noise and reducing the 
influence of environmental cosmic variance on the supernova 
samples. 

Different observations have already produced some 
evidence for shorter time delays: Recent work by 
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Figure 10. Determinations of the SNR in the literature compared with the parametrized versions used in the analysis. The most probable 
t corresponds to the mode of the Bayesian probabilities and the best-fitting model corresponds to the model with the smallest KS test 
D value. The confidence intervals are based on the KS test, corresponding to the models where the non-rejection probabilities a r e 0.32 
and 0.05, for the 68 and 95 per cent confidence intervals, respectivel y. The horizontal err or bars are the bin sizes in lDahlen et al.l i'2004) 
and lBarris fc Tonrvl l2005lV Note that the highest redshift bin from lDahlen et alj 120041) is based on two detected SNe. The best- fitting 
models and confidence intervals are obtained using only the SN sample from S04. 



iMannucci et alJ i2005l) showed that the most efficient host 
galaxies for SNe la production am ong all galaxies are irregu- 
lars, and lDella Valle et al] i2005l) showed that among ellipti- 
cal galaxies it is in particular the radio-loud galaxies, which 
are also bel ieved to be associated with r ecent star forma- 
tion. Also, llGarnavich fc Gallagher! 120051) have shown that 
SNe la of normal luminosity occur particularly in those el- 
liptical galaxies with quite substantial star-formation rates. 
Also, SNe la in galaxy cluste rs seem to indicate time delays 
that are shorter than 2 Gyr jMaoz fc Gal-Yamlliooi) . 

Finally, it is clear that, if the SNe la phenomenon is 
composed of several production channels, all conclusions we 
have drawn apply to the dominant channel. A generally less 
common channel with different characteristics could again 
be the dominant channel in a subset of galaxies with older- 
age or higher-metallicity stellar populations. It is already 
established that under-luminous 91bg-type SNe la prefer- 
ably occur in non-star-forming hosts, such as ellipticals. It 



remains to be explored theoretically, whether the CO WD + 
RG - SD scenario could be related specifically to old popu- 
lations. We should also consider the possibility that we may 
not even have found the dominant prog enitor channel for 
normal SNe la l)Hamuvll2003l: lToudl200a) . 

Even if the SFH peaks at redshift ~ 1 and the recov- 
ered time delays are consequently low, the associated SNR 
(see Fig. 1101 lower panel) tends to be over-estimated in the 
highest-z bin. This effect could be interpreted as the signa- 
ture of a long time delay component that does not contribute 
to the total SNR when the universe is too young. However, 
because the highest-z bin only contains two SNe, this does 
not lead to low non-rejection probabilities. A study of bi- 
modal time delay distributions could only be done with this 
method if the uncertainties in the SFH were significantly 
reduced and the metallicity cutoff on the efficiencies was 
properly quantified. 

If different channels produce SNe la from progeni- 



© 0000 RAS, MNRAS 000, 000-000 



12 F. Forster, C. Wolf, Ph. Podsiadlowski and Z. Han 



tors of distinctly different age or metallicity, then an in- 
crease of the supernova sample could greatly help the 
identification of the various plausible progenitors. How- 
ever, such a data set would most be beneficial if it 
is complemented with host galaxy characterisation (see 
Ivan den Bergh. Li &_ Filippenko| 2005f) and spectra with bet- 



van den Bergh. Li & t UiPDcnko m 
ter signal fsee lBenetti et al] |2005'l. 
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